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ABSTRACT 

We present new 12 CO(J=T0) observations of the barred galaxy NGC 4303 using the Nobeyama 45m 
telescope (NR045) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). 
The Ha images of barred spiral galaxies often show active star formation in spiral arms, but less so 
in bars. We quantify the difference by measuring star formation rate and efficiency at a scale where 
local star formation is spatially resolved. Our CO map covers the central 2f 3 region of the galaxy; the 
combination of NR045 and CARMA provides a high fidelity image, enabling accurate measurements 
of molecular gas surface density. We find that star formation rate and efficiency are twice as high 
in the spiral arms as in the bar. We discuss this difference in the context of the Kennicutt-Schimidt 
(KS) law, which indicates a constant star formation rate at a given gas surface density. The KS law 
breaks down at our native resolution (~ 250 pc), and substantial smoothing (to 500 pc) is necessary 
to reproduce the KS law, although with greater scatter. 

Subject headings: ISM: molecules - galaxies: individual (NGC 4303) - galaxies: ISM - galaxies: spiral 
- radio lines: galaxies 



1. INTRODUCTION 

Barred galaxies often reveal a striking difference in 
star formation activity between the bar and outer sp i- 
ral arms (e.g. IDownes et ai~IH996l: ISheth et al. 112002) . 
Most star forming regions (Ha emission) are associated 
with the outer spiral arms rather than with the bar, even 
when both host a similar amount of molecular gas. 

The difference is often attributed to strong shear mo- 
tions along the bar, by which giant molecular clouds 
(GMCs) could be pulled apart an d disrupted before the 
ignition of star forma tion (e.g. IDownes et al. 1 119961 : 
Schinnerer et al. | 20Q2D . On the other hand, according to 



Koda fc Sofue I (|2006l ). the observed velocity gradient in 



a bar is not large enough to shred gravitationally-bound 
GMCs. 

Unfortunately, previous studies suffered from large un- 
certainties in measurement of gas surface density, due to 
the poor image quality (fidelity) of interferometer obser- 
vations, or to low sensitivity and resolution even when 
the observations from a single dish telescope are included. 
Interferometers have the intrinsic missing flux problem, 
and cannot measure the gas surface density accurately. 
Most interferometric studies covered only a central bar 
region including the nucleus, and did not extend to outer 
spiral arms, because of the limited field-of-view. They 
studies concentrated on nuclear gas c ondensations and 
the ga s fueling galactic centers(e.g. iSakamoto et al. I 
Il999a|) . and did not contrast the difference between the 
bar and outer spiral arms. At low resolution, the veloc- 
ity gradients and gas surface densities appear similar in 
the inner bar and outer spiral arms, perhaps because of 

1 Department of Astronomy, University of Tokyo, Hongo, 
Bunkyo-ku, Tokyo 113-0033, Japan; momo.s.rieko@nao.ac.jp 

2 National Astronomical Observatory of Japan, Mitaka, Tokyo 
181-8588, Japan; sokumura@nro.nao.ac.jp, tsawada@alma.cl 

3 Department of Physics and Astronomy, Stony 
Brook University, Stony Brook, NY 11794-3800, USA; 
Jin.Koda@stonybrook.edu 

4 California Institute of Technology, MS 105-24, Pasadena, 
CA 91125, USA 



beam dilution (e.g. iKuno et al. Il2007l) . 

We quantitatively study the star formation rate (SFR) 
and efficiency (SFE = SFR/Mh 2 ; the ratio of the SFR to 
the molecular gas mass) at a high resolution over the disk 
of the barred galaxy NGC 4303, using new molecular gas 
data at high quality. We capitalize on the recent develop- 
ment of the Combined Array for Research in Millimeter- 
wave Astronomy (CAR MA) interferometer, a nd the 25 
beam receiver (BEARS: ISunada et al. I I2000D and on- 
the-fly (OTF) mapping technique (jSawada et aT~ll2008t l 
at the Nobeyama 45 m radio telescope (NR045). The 
combination of CARMA and NR045 solves the missing 
flux problem, providing the best image fidelity and qual- 
ity, and enabling accurate measurement of gas surface 
density. Our measurement is based on the fundamental 
transition J = 1 — of CO, which ha s often been used to 
deriv e the gas surface density (e.g.. iScoville fc Sanders I 
ll987tlYoung fc Scovillelll991[) . We occasionally refer to 
"molecular gas" as "gas" for simplicity in the rest of the 
paper. 

We also investigate the star formation activ ity in the 
context of the Kennicu tt-Schmidt (KS) law (|Schmidt I 
119591: lKennicuttlll998D . The KS law is the observa- 
tional correlation between the area averaged SFR (Ssfr 
[Mq yr _1 kpc -2 ]) and gas surface density (Sh 2 [Mq 
pc -2 ]) as a power law: Ssfr oc £h 2 - This correlation of 
iKennicuttl (|1998| ) spans a wide range, 5 orders of mag- 
nitude in Sh 2 an d 7 in Esfr, with N = 1.4 ± 0.15 using 
a sample of 61 nearby spiral and 36 starburst galaxies. 

Most recent studies concentrate on the KS law at high 
spatial resolutions, and use local gas and star forma- 
tion rate densities rather than the averages over galac- 
tic disks. The power law index is a topic of interest, 
since it could indicate the physical origin of the cor- 
relation. A linear correlation, N = 1, would indicate 
that the star formation efficiency is a constant over a 
large range of gas surface density; in other words, the 
mere presence of the gas is sufficient for star formation. 
On the other hand, a super- linear correlation, N ^ 1, 
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may suggest that additional conditions or trigger mech- 
anisms, such a s gas dynamics, also play critical roles. 
iWong fc Blitzl (120021) averaged their data azimuthally 
as a function of galactic radius for 7 gas rich galaxies, 
and obtained N = 1.1 and 1.7 ( f or two different extinc- 
tion models). iKennicutt et al. 1 (|2007| ) used the locally 
averaged parameters (520 pc aperture) around Ha and 
24 /mi peaks in M 51, and obtained N = 1.56 ± 0.04. 
iBigiel et all (j2008ft studied the KS law for 18 nearby 
galaxies at 750 pc resolut i on an d derived N = 1.0 ± 0.2. 

Recently. IVerley et all (|2010f ) showed that the KS law 
breaks down at a smaller scale (180 pc). This may in- 
dicate that the mechanism for triggering star formation 
operates at a scale smaller than 500 pc, and that star 
formation is triggered by the environment of GMCs at a 
scale < 500 pc. We study the relation between the gas 
and star formation at a scale of ~ 250 pc. 

NGC 4303 is a nearby barred spiral galaxy with face- 
on geometry and is the ideal target for this study. 
It is a member of th e Virgo cluster (D = 16.1 Mpc; 
iFerrarese et al. I 11996ft , whose bar and spiral arms are 
resolved spatially and kinematically with CARMA (Ta- 
ble [1}. Even a glance at an Ha image is enough to 
recognize the striking difference in star formation ac- 
tivit y between the bar and spiral arms in this galaxy 
(see iKoopmann et al. I [20011 Figure 5). The CO emis- 
sion from NGC 4303 has been observed with various ra- 
dio telesc opes, i.e., the Owens Va lley Radio Observatory 
(OVRO: ISchinnerer et~aT1 [200l. N obey ama Millime- 
ter Array (NMA: lKoda fc Sofue 12006ft. Berke lev-Illinois- 
Maryland A ssociation (BIMA: iHelfer et al. 112003ft , and 
the NR045 dKuno et al. 112007ft . We observe NGC 4303 
with NR045 and CARMA and combine the data to ob- 
tain a high fidelity image. We focus on discussion of star 
formation rate (SFR) and efficiency (SFE) as a tracer of 
star formation activity at high spatial resolution within 
the galaxy. 

The outline of this paper is as follows. In §2, we present 
our CO(J=T0) observations and data analysis. Results 
are described in §3. SFR and SFE in the bar and spiral 
arms of NGC 4303 are discussed in §4. The KS law is 
discussed in §5. A summary is given in §6. 

2. OBSERVATIONS 

We observed the 12 CO(J=1-0) line emission (/ rcs t = 
115.2712 GHz) toward the barred spiral galaxy NGC 
4303 with NR045 and CARMA. The combination of 
NR045 and CARMA interferometer data provides an 
unprecedented level of image fidelity at a high resolu- 
tion. Interferometers have an intrinsic limitation: they 
cannot obtain the shortest spacing components in u — v 
plane, and thus are not sensitive to the most extended 
structures. The single-dish NR045 covers the short spac- 
ing data. The combination of the two telescopes is ideal 
for the study of extended nearby galaxies, enabling accu- 
rate measurements of surface brightnesses and gas sur- 
face densities. The observing parameters are listed in 
Table 2. 

2.1. NR045 

The single dish NR045 observations were carried out in 
March 2008 at the Nobeyama Ra dio Observatory. Th e 
twenty five-beam array, BEARS (jSunada et al. I [2000h . 



consists of double-side band (DSB) superconductor- 
insulator-supcrconductor (SIS) receivers on a 5 x 5 grid 
with a separation of ~ 41'.'1 at 115 GHz. We used the 
digital auto-correlator for the NR045 (jSorai et al. 112000ft 
with a bandwidth of 512 MHz (~ 1340 km s" 1 ) and 500 
kHz resolution (1.3 km s _1 at 115 GHz). 

The chopper wheel method was used for atmospheric 
and antenna ohmic loss corrections and to calibrate the 
antenna temperature in the double-sideband T^(DSB). 
A typical system temperature T sya was 300 - 400 K 
in the DSB during our observations. T^(DSB) was 
converted to the single sideband (SSB) temperature 
T^(SSB) by applying scale factors (hereafter, we use T* A 
to refer to T^(SSB)). The scale factors were measured 
by comparing images of the Orion IRc2 region in the 
CO(J=1-0) line obtained using BEARS and the single- 
beam SSB receiver S100. The main beam efficiency of 
NR045 was ?7mb = 0.30 at 115 GHz. The main beam 
temperature is calculated as Tmb = T a /t]mb- 

The telescope pointing was checked roughly every 30 
minutes by observing the nearby SiO maser source R-Crv 
with another SIS receiver, S40 near 40 GHz. The stan- 
dard deviation of the pointing offset between the optical 
axes of S40 and BEARS was 0'.'87, which is measured 
by the observatory. The pointing error was less than 6" 
during our observations. We mapped the 8' x 8' region 
centered on NGC 43 03 with the OTF observing mode 
(jSawada et al. 112008ft . which corresponds to a 37 kpc x 
37 kpc area at the distance of NGC 4303. This area cov- 
ers the optical extent of the galaxy. We alternated be- 
tween continuous scans in RA and DEC, as systematic 
noise appears dominantly along the direction perpendic- 
ular to the scan, which causes it to cancel out when the 
scans are combined. The 25 beams draw a regular 5" 
stripe in a single scan, which is finer than the Nyquist 
spacing of the 15" NR045 beam — no spatial informa- 
tion is missed, which is critical in producing good u — v 
coverage and in combining with CARMA data. We con- 
struct data cubes only for the 5' x 5' area covered by all 
25 beams (Figure [1]). 

We used the data reduction package NOSTAR 
(Nobeyama OTF Software Tools for Analysis a nd Re- 
duction) developed at NRO (jSawada et al. (12008ft . After 
the DSB-to-SSB conversion, we subtracted spectral base- 
lines with a first order polynomial fit and flagged bad in- 
tegrations. The RA and DEC scans provide data points 
roughly along a 5" lattice. We used a spheroidal func- 
tion to re-grid them onto a regular grid to produce an 
integrated data cube. The Spheroidal function has low 
sidelobes in the Fourier domain and is advantageous in 
de-convolving the map to produce u — v data. The con- 
volution smears out the resolution of the map produced 
by the 15" NR045 beam to 22'.' 1. We produced separate 
data cubes for the scans taken in the RA and DEC direc- 
tions and combine them using the basket-weave method 
in NOSTAR. Lower weights were applied to the data at 
small wave-numbers in the scan directions in the Fourier 
domain, which significantly reduces the systematic noise 
in the combined map. 

The final data cube has a spatial resolution of 22" 1 and 
a pixel scale of 7"5, which is the half size of 15" NR045 
beam. The rms noise is 75 mK at the velocity resolution 
of AV = 2.6 km s" 1 . 
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2.2. CARMA 

We observed the central 2'3 region of NGC 4303 with 
CARMA in the C-configuration (October and November 
2007) and in the D-configuration (June 2008). CARMA 
is a 15-element array, which consists of nine 6.1-m anten- 
nas (formerly BIMA) and six 10.4-m antennas (formerly 
OVRO). The angular resolution and baseline coverage of 
C-configuration (D-configuration) is 2" (5") and 26 - 370 
m (11 - 148 m). The spatial resolution corresponding to 
the minimum baseline (11 m) is 49". We adopted a 19 
point hexagonal mosaic pattern with 30" spacing (i.e., 
the Nyquist spacing of a 10.4-m primary beam; Figure 
[I]). We used a 3 mm SIS receivers, and digital correlators 
with 3 bands, each of which has 62 MHz of bandwidth in 
a single side band with 977 kHz resolution (2.5 km s -1 
at 115 GHz). 

The bright quasar, 3C273, was observed for bandpass 
and gain calibrations. We observed it roughly every 20 
minutes to calibrate for atmospheric and instrumental 
gain variations. We corrected the telescope pointing by 
observing 3C273 with the 3 mm receiver. We frequently 
monitored and corrected the pointing during tr acks, us- 
ing th e optical pointing method developed by iCorder I 
(|2008l ) , which involves observing bright optical stars close 
to 3C273 every gain calibration cycle (around twenty 
minutes). The pointing offset vector between the radio 
and optical pointings was measured at the beginning of 
the observations. 

We used the data reduction package MIRIAD 
(|Sault et al. 1 119951 ) for calibration and data analysis. 
The rms noise level of the CARMA data is 43 mJy 
beam -1 in a channel of Ay = 2.6 km s" 1 . The syn- 
thesized beam size is 3"1 x 2"5. 

2.3. Combining NR045 and CARMA data in u — v 

space 

We combined CA RMA and NRQ45 da t a in the Fourier 
domai n as done in iKoda et al. I (|2009f) ; iKurono et al. I 
([2009(1 . The NR045 provided the short spacing data 
which filled the central hole in the CARMA u — v cover- 
age. The combined image achieved higher angular resolu- 
tion than the NR045 map and recovered the flux missing 
in the CARMA map. 

By utilizing MIRIAD^ and new routines that we 
developed, we combined NR045 and CARMA data in 
the following way: 

1. Deconvolution of single dish map with single-dish 
beam 

The point spread function (PSF) of the NR045 
map is a convolution of the NR045 beam pattern 
and the spheroidal function that we applied to re- 
grid the data points and create the final data cube. 
The map should be deconvolved by the PSF be- 
fore being Fourier-transformed to the u — v plane. 
The NR045 beam size is Gaussian with an effective 
beam size of 22" 1. 

2. Multiplying dummy 2' Gaussian primary beam 
This step produced the same data structure for the 
single dish data as is used by the interferometer 



data in order for the files to be read into MIRIAD. 
We used a 19 point hexagonal mosaic technique 
in our CARMA observations, so the NR045 map 
was divided into 19 smaller maps corresponding 
to CARMA's mosaic points. The pseudo-primary 
beams were multiplied by each of the 19 maps to 
mimic CARMA observations. CARMA has three 
types of primary beam pattern with a FWHM of 
64" for 10.4 m - 10.4 m, 85" for 6.1 m - 10.1 
m, 115" for 6.1m - 6.1 m at 115 GHz. To mimic 
this behavior, we applied a Gaussian primary beam 
with a FWHM of 2' to the deconvolved NR045 
map at each of the 19 pointings and produced 19 
maps. The Fourier transformation of these maps 
are the u — v data we need (see step 4) . 

3. Determining weightings 

The relative weighting of NR045 and CARMA 
data changes the synthesized beam pattern. By 
changing the weighting, we can control the side- 
lobes of the synthesized beam and the noise 
level. This is similar to the conventional weighting 
schemes (e.g. natural and uniform weighting, or ro- 
bust parameter). We decided to use the weighting 
that minimized the sidclobes. 

4. Constructing visibility data from NR045 data 
All of the visibility data of NR045 is built after ap- 
plying the relative weight balance b etween the visi- 
bility data of NR045 and CARMA (jKurono et al. I 
I2009D . We made 57 visibility data sets, i.e., for 19 
pointings x 3 primary beams. All 57 NR045 vis- 
ibility data sets are concatenated to produce the 
final NR045 u — v data. 

5. Imaging 

We obtain the combined data cube by adding to- 
gether the visibility data sets from NR045 and 
CARMA, and by Fourier transforming the result, 
we created a combined image. The total flux of the 
combined image is similar to that of the NR045 im- 
age. The rms noise of the combined cube is 34 mJy 
beam -1 with AV = 2.6 km s . The synthesized 
beam is 3"2 x 3"2. 

Bright CO emission is detected throughout the disk of 
NGC 4303, even in the inter-arm region in NR045 map 
(Figure H (a)). On the other hand, the CARMA map 
shows emission only from a compact component in the 
nuclear region and the spiral arms (Figure [2] (b)). The 
combined data map clearly shows strong CO emission 
in the nucleus, the bar, and the disk region along the 
two arms (Figure [5] (c)). A comparison of the total flux 
in the NR045, CARMA, and combined images clearly 
shows that the missing flux has been completely recov- 
ered (Table EJ). 

2.4. Archival data 

We obtained archival Ha dKnapen et al~ll2004[ ) and K- 
band ([Mollenhoff fc Heidt Il200lf ) images from the NASA 
Extragalactic Database (NEE0) and a 24 /zm image (P.I. 
Kotak) from the Spitzer Space Telescope data archive 
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using the Leopard software. The retrieved 24 (im image 
suffered from residual sky background emission, which 
we subtracted using the Image Reduction and Analysis 
Facility (7iL4i0) package. 

UV photons from recently-formed, young, massive 
stars primarily ionize the surrounding gas (H II regions) 
or heat dust grains (thermal radiation). Thus, the combi- 
nation of Ha and 24 /im dust continuum emissions traces 
the s t ar formation rate fairly accurately ([Calzetti et al. I 
120071 iKennicutt et al. 1 12007ft. The SFR i s calcu lated 
using the formulae f rom iKennicutt et al. I (f2007l ) and 
ICalzetti et aTl (120071 ). namely 

SFR[M year" 1 ] = 7.9 x 10~ 42 L(Ha) corr [crgs s" 1 ] 



7.9 x 10~ 42 (L(Ha) ob 



(0.031 ± 0.006)L(24^m))[erg s" 

(1) 



where L(Ha) and L(24 /im) are the Ha and 24 /im lumi- 
nosities. From the Ha and 24 fim images, we found that 
Ha emission defined the sites of star forming regions very 
well, although it does not quantitatively trace the star 
formation rate (SFR). Since the Ha image has a much 
higher resolution (1'.'5) than the 24 fan image, we use the 
Ha image in discussing the offsets between star forming 
regions and parental gas qualitatively. For quantitative 
measurements, we convolve the Ha data, and smooth the 
Ha and CO images to match the spatial resolution of the 
24 /im image (6" - 500 pc). We clip the SFR and CO 
data at the 2a r.m.s. noise level. The SFR map is shown 
in Figure E] (left). 

3. RESULTS 

3.1. Molecular Gas Distribution 

Molecular gas distributions in barred spiral galaxies 
have been studied by many authors, fro m high resolu- 
tion observations of central bars (e.g. iKennev et al. I 
ll993HIshizukTlll999HSakamoto et al. lll99 9a.b) t o lower 
resolu tion observations of larger disks (jSheth et al. I 
120001 ) . The gas distribution in NGC 4303 is typical 
among barred spiral galaxies; it shows two narrow ridges 
(namely, offset ridges) at the leading side of the bar, and 
spiral arms in the outer part, and coincides with dust 
lanes in optical images. Figure [5] (d) shows an overlay of 
CO contours on an archival HST F450 band image (P.I. 
Smartt). 

The CO emission is concentrated in the central 40" (~ 
3.2 kpc), which includes the offset ridges and circumnu- 
clear disk. In particular, the circumnuclear disk shows 
a circular structure (r = 4'.' 4 ~ 350 pc), extending from 
1541.3 km s~ 4 southeast of the nucleus to 1611.5 km s _1 
to the northwest in the velocity channel maps (Figure 0] 
FigureEJ . This structure is calle d a ring or part of the off- 
set ridges (i.e., inner s piral arms; iSchinnerer et al. 12002b 
iKoda fc Sofue I [20061 ) . The circumnuclea r disk hosts a 
spira l arm in the UV continuum image (|Colina et al. I 
119971) . and the CO ring/spiral arm structure is offset from 
the UV spiral arm. 

Low-level CO emission is also present in the space out- 
side of the offset ridges (namely, inter-offset ridge re- 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



gions). Strong concentrations of gas exist at the outer 
ends of the offset ridges (Figure [2] (c)). Molecular spiral 
arms extend from those end-points toward the outer re- 
gions. The map resolves the widths of the offset ridges 
of the bar and the outer spiral arms at high resolution. 
The surface densities of the gas in the offset ridges and 
outer spiral arms are not similar at the high resolution, 
differing by factor of ~ 2. 

We derive the inclination and position angle of the 
galaxy using the task GAL in the AIPS data reduction 
package. The results are shown in Table 1. We use these 
values for the rest of our analysis. 

3.2. Molecular Gas Mass and Surface Density 

], We estimate the gas mass from the CO integrated in- 
tensity Sco ( see Figure 0(c)), using the Galactic J(CO)- 
to-V(H2) conversion factor Xco- We adopt the standard 
value X C o = 2.0 x 10 20 cm" 2 [K km s" 1 ]" 1 (jStrong I 
[19871) . The gas mass is then 



M H , =1-0 x 10 4 



s, 



CO 



Jy km s 1 



(— V 

VMpc J 



3.0 x 10 20 cm~ 2 (K km s" 



(2) 



©i 



where D is th e distance to the galax y , which assumed to 
be 16.1 Mpc (jFerrarese et al. Ill996| ). The total gas mass 
of the galaxy is calculated from the NR045 map to be 
5.3 x 10 9 M Q . The CARMA mosaic covered the circular 
area with a radius of 80". The enclosed mass is 3.4 x 10 9 
M Q . The average gas surface density in the central 160" 
of the molecular disk is 36 M Q pc~ 2 . 

4. MOLECULAR GAS AND STAR FORMATION 
4.1. Spatially Resolved Sites of Star Formation 

Both CO and Ha emission appear clearly along stel- 
lar spiral arms, and the Ha arms are offset toward the 
leading side of the CO arms. Similar offsets are evident 
between CO and 24 /im emission. The 24 fim emission is 
at the leading side of the CO emission. These offsets are 
due to the delay of star for mation from the acc umulation 
of gas in the spiral arms (|Egusa et al. Il200l ; approxi- 
mately, an offset would be (fi — Q, p )tsFi where and 
fi p are the angular speed of the gas and pattern speed 
of spiral arms, respectively, and t$F is the timescale of 
star formation. The sites of star formation are resolved 
in our observations. The analysis of local star formation 
rate and efficiency should take this into account, since 
the parental gas and star forming regions would not be 
correlated at the high spatial resolution. This means that 
the KS law breaks down at this resolution. We need to 
smooth the data before characterizing the spatial varia- 
tion of the star formation efficiency in §4.2. 

Figure [S] (a) , (b) show a close-up of the offset between 
CO (black contours) and Ha emissions (greyscale). In 
the upper panels, the arrows indicate the direction of gas 
flow, assuming pure azimuthal rotation and a time lag 
between star formation and gas accumulation in spiral 
arms. Figure [6] (c) shows the phase diagrams of CO 
(contours) and Ha (greyscale). Again, the Ha emission 
is found predominantly at the leading side of the CO 
spiral arms over the entire range of radii (over 40", ~ 3.1 
kpc). 
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The offsets between the CO and Ha arms range from 
to 1 kpc, with an average of ~ 500 pc at radius, r > 
40" (~ 3.1 kpc; i.e., the area we define as outer spiral 
arms in §4.4). Thus, we will smooth the data to the 
spatial resolution of ~ 500 pc for later discussions; it also 
especially affects the analysis of SFE. We will discuss this 
f urther in § 4.4 . 

lEgusa et al. I (|2004[ ) introduced a new method to es- 
timate the timescale of star formation using these off- 
sets, measuring the offset /tim e lag between ga s accu - 
mulation and star formation. iTamburro et al. 1 (|2008T ) 
applied this method to H I and 24 (im data. Their off- 
sets < 250 pc (calculated from their paper), are trans- 
lated to the t i mesca le of ~ 1 Myr. On the other hand, 
lEgusa et al. 1 ((20091) used CO and Ha data for NGC 
4303, and derived the offsets of 700-1000 pc and the 
timescal e of ~ 10 Myr, an ord er of magnitude longer 
than the ITamburro et all (120081 ) resul ts. Our offsets are 
closer to those in lEgusa et al. I (|2009f ) . To understand 
this discrepancy, we compared 24 /xm. and Ha images, 
and found no significant offsets between the two. Thus, 
Ha traces the locations of star forming regions as well 
as 24 /im. The discrepancy may come from the differ- 
ence of CO and H I. If the H I emission traces the gas 
dissociated by star formation, and not the natal gas for 
star formation, it explains the small offsets between H I 
and 24 /im emissions. Our simple analysis with CO and 
Ha/24 (jm favors th e longer star formation timescale of 
lEgusa et all ([20091) in case of NGC 4303. These results 
should be confirmed with a larger sample of galaxies. 

4.2. Star Formation Rate and Efficiency 

We compare the star formation rate (SFR) and star 
formation efficiency (SFE) between the bar and outer 
spiral arms, as well as within other components (i.e., cir- 
cumnuclear disk, inter- arm regions), at the spatial reso- 
lution of - 500 pc. The SFE is the ratio of the SFR to 
the gas mass within an aperture, namely 



SFElyr- 1 ] = 



SFR[M Q yr- 
M Ha [M ] 



(3) 



where Mh 2 is the molecular hydrogen gas mass. We cal- 
culate the SFR and gas mass within an aperture and use 
the averaged parameters to derive the SFE. We smooth 
our CO and Ha images to the same resolution as the 
24 /xm image (6" resolution) in order to correlate the 
parental gas and star forming regions. This resolution 
corresponds to a linear siz e of ~ 500 pc, similar t o the 
linear resolution at which iKennicutt et al. 1 (|2007l) find 
the correlation between gas surface density and associ- 
ated star formation rate. 

A SFR map is shown in Figure [3] (a). High SFRs 
are apparent in intense star forming regions which are 
sparsely distributed along the spiral arms, as well as 
in the circumnuclear disk. However, SFRs are lower in 
the bar, even though the gas surface densities along the 
bar offset ridges are as high as those in the spiral arms. 
The low SFR in the bar has been suggested qualitatively 
by pr evious studies (most systematically bv lSheth et al. I 
120021) , and we discuss this quantitatively in the following 
sections. 

Figure [3] (b) shows a map of the SFE. Low SFEs are 
seen in the bar (especially around the offset ridges) and in 



the inter-arm regions. Exceptionally high SFEs appear 
at the leading side of the CO arms. This indicates that 
the offsets between Ha and CO emission are not fully 
smoothed out even at the resolution of ~ 500 pc (see § 
3.2), even t hough the parental gas s hows a fair correlation 
with SFR (|Kenriicutt et al. 11200% . 

4.3. Azimuthally Averaged Star Formation Efficiency 

We first compare the gas surface density (Eh 2 ), area 
averaged SFR (£sfr), and SFE between the bar and spi- 
ral arms on the basis of azimuthal averages (Figure [7J . 
The average values are derived in 1 kpc (12'.'8) width an- 
nuli, which is enough to average out the spatial offsets 
(see § 3.2). The lkpc resolution can still separate the 
bar, spiral arms, and nucleus without significant contam- 
ination from each other. Figure [3] (a) shows the defini- 
tions of nucleus, bar, and outer spiral arms. The nucleus 
might include a portion of the bar region; however, we 
are mainly interested in the differences between the bar 
and spiral arms, so it does not affect our analysis. We 
will comment on the properties of the circumnuclear star 
formation in § 4.4. 

Figure [7] shows the radial distributions of the area aver- 
aged SFR Ssfr (a), gas surface density Sh 2 (b), and SFE 
(c). The definitions of the nucleus, bar, and spiral arms 
are indicated with arrows in Figure [Jj (a) . SFE drops 
abruptly from the circumnuclear disk to the bar, and 
increases toward the spiral arms (Figure [7J (c)). Quan- 
titatively, the average SFE is about twice as high in the 
spiral arms as in the bar. 

4.4. Locally Averaged Star Formation Efficiency 

The star formation efficiency (SFE, Eq. [3|) is a correla- 
tion between star forming regions and their parental gas. 
We have resolved the scale where the star forming re- 
gions and gas are spatially separated; therefore, we have 
to smo oth the data to derive a SFE with physical impor- 
tance. IKennicutt et al. 1 (|2007t) found a good correlation 
between the surface densities of SFR and gas mass at the 
500 pc scale. This is similar to the typical separation of 
Ha and CO emission that we see in NGC 4303 (EH- 
We smooth Ha and our CO images to the scale of the 
Spitzer 24 fj,m image (6"; ~ 500 pc), and compare local 
SFR and SFE between the bar and spiral arms at this 
scale. 

Four regions (i.e., bar, spiral arms, inter-arms regions, 
and circumnuclear region) are defined as in Figure [8] 
(b). The circumnuclear region is a 6" diameter circle 
(~ 500 pc), based on a near- infrared (NIR) image of 
Perez-Ramfrez et al. (2000). We define the bar as an 
ellipse with the deprojected major and minor axis diam- 
eters of 40" x 26" (3.1 x 2.0 kpc), based on a K-band 
image. This de finition yields an ell ipticity similar to the 
one derived by lLaine et al. I (|2002l ) based on a near in- 
frared ii-band image, and a m ajor axis diameter consis- 
tent with the determination of iSchinnerer et al. I (|2002|) . 
The spiral arms and inter-arm regions are separated by 
eye using the K-band image. Figure [8] (b) shows our 
definitions overlaid on the i4f-band image. 

Figurc[5]shows the histograms of Sh 2 (left), Ssfr (cen- 
ter), and SFE (right) for the spiral arms and bar (top and 
bottom). High gas surface density regions are present 
only in the bar, but not in the spiral arms. However, high 
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SFR and SFE regions appear only in the spiral arms and 
are absent in the bar. Therefore, high SFR and SFE are 
evident in the spiral arms in the local average, as well as 
in the azimuthal average. 

We calculate the excesses of SFR and SFE in each re- 
gion with respect to the averages over the entire disk 
excluding the nucleus (Table 4). The SFR and SFE in 
the bar are 10 % and 30 % smaller than the disk aver- 
ages, respectively. Star formation activity in the bar is 
low compared to the entire disk. On the contrary, the 
SFR and SFE in the spiral arms are about 10 % and 40 
% higher than those over the disk, respectively. Active 
star formation preferentially occurs in the spiral arms. 
SFR and SFE are about 30 % and twice as high in the 
spiral arms as in the bar. The results confirm quantita- 
tively the notion that star formation in bar is less than 
in spiral arms. 

The effect of the spatial offsets between Ha and CO 
emissions is significantly reduced by the smoothing to es- 
timate local SFE (~ 500 pc resolution). However, it still 
shows some residual effects in SFEs, since the average of 
this offset is also ~ 500 pc. If the offset is below 500 
pc, parent GMCs and star forming regions exist within 
our resolution for SFE discussions. If the offset is above 
500 pc, we cannot recognize both parent GMCs and star 
forming regions within the resolution. SFEs obtained 
in such a region will not fairly represent the local SFE. 
This produces SFEs above 10 -78 yr _1 are the artifacts 
(Figure [3]). 

4.5. Circumnuclear Star formation 

The SFR and SFE are the highest in the nuclear region. 
The azimuthal averages of Esfr and SFE (Figure H (a), 
(c)) clearly show this tendency, being about 5 times and 
twice higher than the disk averages, respectively. The 
locally averaged SFE in the nucleus is also high. The 
locally averaged SFR and SFE arc higher by a factor of 
6.9, and 1.3, respectively than those of the entire disk 
region. 

We, however, note that NGC 4303 is suggested 
to harb or an AGN in addi ti on to a starburst 
(e.g. IKennicutt et aT~l 119891: iColina et al~l 119971 : 
iTschoke et al. II2000D . and our analysis could suffer from 
AGN contamination. We checked the spectral energy 
distribution (SED) of the nucleus in 1-24 /im, using 
the archival data ( Catalog of Infrared Observation^, 
iGezari et al. 111993] ) . We found a flat spectrum from 10- 
24 /im , similar to AGN spectra ([Alonso-Herrero et al. I 
(2003])). Thus, 24 /im emission from the nucleus is dom- 
inated by the AGN. This indicates that the much higher 
values of SFR in the nucleus are overestimated. 

4.6. Comments on metallicity- dependent Xco 

We adopt the standard conversion factor Xco- 
However, there are some arg uments supporting a 
metallicity-depende nt Xco (e.g. lArimoto et al. I [l996t 
iBoselli et al~l l2002h . NG C 4303 has a relative ly large 
radial metallicity gradient (jSkillman et al. Ill996|) , which 
may affect the result that SFE is low in the bar (inner 
disk) and high in spiral arms (outer disk) . We confirmed 
as follows that this does not change the results (except 
for the nucleus, where metallicity is exceptionally high). 

8 http:/ /ircatalog. gsfc.nasa.gov/ 



The average metallicities (12 + log[0/H]) in the nu- 
cleus, bar, and spiral arms (the defini tion is same as 
Figure M (b)) are 9.5, 9.3, and 9.2 from iSkillman et al. I 
(11996ft. If we ado pt the Aco-metallicity relation of 
IBoselli et al. I ()2002l ). the gas surface densities decrease, 
and the average SFEs increase (by a factor of 4.2, 2.3, 
and 1.8) respectively. The difference in SFE between the 
bar and spiral arms is reduced only by 30 %. This is too 
small to change the trend that the SFE is roughly a fac- 
tor of two higher in spiral arms than in bar. A significant 
effect appears only in the nucleus, where a much higher 
metallicity results in a significantly higher SFE (a factor 
of 2 - 4 higher than that from the standard Aco)- 

5. THE KENNICUTT-SCHMIDT LAW 

The Kennicut-Schmidt (KS) law shows a tight correla- 
tion betwe en the gas surfac e density E gas and SFR den- 
sity Esfr (jKennicutt 1119981 ) as we mentioned in § 1. The 
difference in SFE between the bar and spiral arms may 
affect this correlation. Figure [10] (a) shows the surface 
densities of gas (Eh 2 [M© pc -2 ]) and area averaged star 
formation rate (Esfr [M© yr -1 pc -2 ]) for our data at ~ 
500 pc resoluti on. We use t he sam e axis ranges as those 
in Figure 9 of IKennicutt I (|1998l ). spanning five orders 
of magnitude in Ejj 2 and seven orders of magnitude in 
Esfr, a range which covers from normal spiral galaxies 
to circumstellar starburst. 

The central region of our data could be contaminated 
with AGN, and the SFR calculated form 24 /im could 
be overestimated. Therefore, we fit straight lines to two 
sets of data, i.e., with and without the central region, 
to see the effect of AGN on the KS fits. Figure [TD] (a) 
shows our least square fit results: the orange and blue 
solid lines are for the data with and without the cen- 
tral region, and are Esfr = 10- 2 - 61±004 x vj^is±o.02 

and E SF R = lO" 178 ^ 05 x £°^ 7±0 - 03 , respectively. The 
contamination of A GN appare n tly in creases the index 
N. The analysis of IKennicutt I (|1998|) might have been 
contaminated with AGNs, which are typically found in 
galaxies with circumnuclear starbursts. 

5.1. Comparison with Previous Studies 

Figure [TU] (a) compares our result with those from 
previous studies. The solid black line is the result of 
IKennicutt I (|1998l ): N = 1.4 ± 0.15. They calculated gas 
surface densities with CO data (gas density), and H I and 
CO data (total gas density). The grey region indicates 
the scatter in Ke nnicutt (1998, their Figu re 9). Dashed 
lines are the fit in IKennicutt et a"l~l (|2007l ): the large and 
small dashes are the fits for H I and CO (N = 1.56 ± 
0.04), and for CO-only (N = 1.37 ± 03) respectively. 
The dotted line is from iBigiel et al. I (|2008| ). and is for 
N = 0.96 ± 0.07 Our res ult, N = 6 7, is closer to 
IBigiel et all (|2008t) . than to iKennicuttl (|1998l) . though 
it's smaller than the linear correlation N =1. 

Our data lie roughly in the range of Kennicutt (1998, 
grey region), although some data points deviate from the 
exact range. The deviations occur on the spiral arms, 
where the offsets between the gas and star forming re- 
gions are large (§4.4). At the 500 pc resolution, the KS 
law is present, but suffers slightly from the offset effect. 
Two types of scatter appear in Figure [10] (a), (b). One 
is the scatter of NGC 4303 itself, i.e., its deviation from 
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the average of galaxies, and the other is within the disk 
of NGC 4303. The distribution of our points is concen- 
trated in a small regio n, and the scatter w ithin the galaxy 
is smaller than that of iKennicutt I ()1998f ). Therefore, the 
scatter among galaxi es seems dominan t. Overall, NGC 
4303 fits on the plot of iKennicuttl (| 19981 ). but has slightly 
higher SFE than the average. 

Both Eh 2 and Ssfr show one order of magnitude scat- 
ters in Figure [10] (b), which is consistent with the results 
of Bigiel et al. (2008, their Figure 4). Therefore, an 
order of magnitude scatter of star forma t ion ac tivity is 
common at a sub-kpc scale. ILerov et al. I (|2008D showed 
that SFE varies strongly with local conditions. Since the 
KS law breaks down at the scale of a few 100 pc, we could 
attribute the variation of SFE to local environments at 
the scale of a few 100 pc. One of the key environmental 
factors might be galactic dynamics (e.g. spiral arms). 

We verify a scatter of NGC 4303 itself among nearby 
galaxies, i.e. a scatter from the KS law. Rough estima- 
tion of scatter shows that our results are a factor of ~ 
5 higher than the that of IKennicutt I ([1998D . However, 
the KS plot has a ± 1 order of magnitude scatter, as we 
discussed in §5.1. Even though our data are the fit by 
IKennicutt"! (|1998D , they are still within the scatter (grey 
region). 

6. SUMMARY 

We observed the barred spiral galaxy NGC 4303 
in the 12 CO(J=1-0) line with NR045 and CARMA. 
The combination of NR045 and CARMA provided an 
unprecedented high image fidelity as well as a high 
angular resolution (3'.'2 ~ 250 pc), which arc critical for 
the accurate measurements of gas surface density and 
mass at high resolution. We discussed SFR and SFE 
quantitatively. Our results are summarized as follows: 

1. CO emission is detected over the entire disk, i.e., 
almost everywhere including interarm regions and 
the downstream side of the bar. There are remark- 
able concentrations along the offset ridges of the 
bar and in the ring structure in the nucleus (r ~ 1.6 
kpc area). The gas in the spiral arms extend from 
the end of the offset ridges toward the outer region. 
The surface densities in the outer spiral arms and 
offset ridges are similar at high resolution. 

2. Spatial offsets between Ha and CO peaks exist 
along the spiral arms. Ha emission is seen at the 
downstream side of gas flow, while the CO emission 
is upstream of the gas flow. The delay of star for- 
mation from the formation of GMC on spiral arms 
would cause such offsets. 

3. The azimuthal averaged SFE decreases steeply 
from the circumnuclear disk to the bar, and in- 
creases toward the spiral arms. The comparison of 



SFE in the bar and spiral arms shows that SFE is 
about twice as high in the arms as those of in the 
bar. 

4. Extreme Esfr and SFE are found in the spiral 
arms, but not in the bar, indicating that the trigger 
of star formation is related not only to the amount 
of available gas, but also to the environment, such 
as galactic dynamics around spiral arms and the 
bar. The presence of the active star forming re- 
gions along the spiral arms confirms the visual im- 
pression that star formation is more active in spiral 
arms, or reduced significantly in bar. 

5. The SFE derived with a metallicity-dependent 
Xqo does not change the conclusion, i.e. higher 
SFE in the spiral arms than in the bar, since the 
difference between the bar and the spiral arms is 
reduced by only around 30%. However, SFE in the 
circumnuclear regions is a factor of 2 - 3 higher 
than the results with the SFE derived by a stan- 
dard Xqo ■ since mctallicity of the circumnuclear is 
significantly high. 

6. The KS law appears to break down at our highest 
spatial resolution (~ 250 pc); due to the spatial off- 
sets we find between parental gas and star forming 
regions. The correlation reappears if we smooth the 
images to ~ 500 pc resolution. The least squares fit 

of all data in the KS law is S SFR = l(r 2 - 61 ± 004 x 
s i.i8±o.02_ NGC 4 3Q3 lieg roughly in the range 

of IKennicutt I (|1998fl . though it has slightly higher 
SFE than average value among nearby galaxies. 
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Fig. 1. — : Optical image of NGC 4303 overlaid with regions observed using the NR045 (blue outline) and CARMA 
(red circles). The 25 BEARS beams are spread evenly inside the region outlined in green, which we used for data 
reduction. 
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Fig. 2.—: Integrated intensity map for (a) NR045, (b) CARMA, (c) NR045 + CARMA (combined data), and (d) 
combined data overlaid on HST F450 image. The contours trace 10% intervals from 10% to 90% of the peak, la noise 
levels are 1.5 K km s -1 (NR045), 810 mJy beam -1 km s" 1 (CARMA), and 710 mJy beam -1 km s -1 (combined data). 
The peak intensities of the NR045, CARMA and NR045 + CARMA images arc 15.7 K km s -1 , 13.0 Jy beam -1 km 
s _1 , and 20.8 Jy beam -1 km s _1 , respectively. The NR045 + CARMA image achieves a high resolution as well as 
recovering extended emission components. The CO emission overlaps with dust lanes, which are seen as extinction in 
the HST image. Red circles arc the central 2 '3 region observed with CARMA. 
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Fig. 4. — : Channel maps of the combined data from 1475.0 km s 1 to 1566.0 km s 1 . The contours are 3, 6, 9a, 
where the ler noise level is 34 mJy beam^ 1 . 



TABLE 1: Parameters of NGC 4303. 



Parameters 



Value 



Name 

Hubble typc a a 

v ist .ib 

Distance c c 

Linear scale of l" 

Diameter d a 

Major 

Minor 

Position Angle d 
Inclination f d 
AGN type" 



NGC 4303 / M 61 
SABbc 

1556.5 km s" 1 
16.1 Mpc 
78.1 pc 

6.'5 
5.'8 
312.2° 
27.8° 

Sy2 / H II 



de Vaucouleurs et al. (1991) Koda & Sofue (2006) Ferrarese et al. (1996) This study. 



b a 
h b 



hi 
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TABLE 2: NRQ45 and CARMA observational parameters. 



Parameters 


NR045 


CARMA 


Date 


March 9th to 12th, 2008 


October to November, 2007, June, 2008 


Field center (J2000) 


12 ,l 21 m 54.90 s , +4 d 28 m 25.1 s 


12 h 21 m 54.9 s , +4 d 28 m 25.6 s 


Observed Frequency 


114.67641 GHz 


114.67641 GHz 


Observed Region 


8' x 8' (37 kpc x 37 kpc) 


19 points mosaicing 


Bandwidth 


512 MHz 


3 X 62 MHz 



TABLE 3: NR045, CARMA and combined data results. 



Parameters 


NR045 


CARMA 


Combined data 


Velocity resolution AV 


2.6 km s" 1 


2.6 km s" 1 


2.6 km s -i 


The rms noise in a channel 


73 mK 


43 mjy beam 


34 mjy beam -1 


The rms noise in an intensity map 


1.5 K km s- 1 


810 mjy beam -1 km s _1 


710 mjy beam" 1 km s — 1 


Beam / Synthesized beam size 


22'.' 1 


3'.'1 X 2'.' 5 


3'.' 2 X 3'.' 2 


Total flux 2.£ 


X 10 3 Jy km s" 1 


9.8 X 10 2 Jy km s" 1 


2.9 X 10 3 Jy km s" 1 



TABLE 4: The average of SFR and SFE over the entire disk regions, bar and spiral arms. 





The entire disk 


The bar 


The spiral 


arms 


SFR [M yr-^j 


8.3 x lO - * 


7.6 x 10" 2 


9.8 X 


10~ 2 


SFE [yr- 1 ] 


5.7 x 10~ 9 


4.2 x 10~ 9 


7.4 x 


10 -9 
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Fig. 6. — : (a), (b) The distributions of CO (black contours) and Ha (greyscalc) emission in two regions of NGC 4303 
are shown. The direction of gas flow is determined assuming a circular galactic rotation and that the young massive 
stars (traced by Ha) flow downstream from the molecular spiral arms where they were formed, (c) The phase diagram 
of CO (contours) and Ha (grey scale). 
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Fig. 7. — : Radial distributions of (a) Ssfr; (b) £h 2 ! ( c ) SFE. The error bars are one r.m.s. 
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Fig. 8. — : Definitions of the nucleus, bar, and spiral arm (outer disk) regions. Left: azimuthal definitions. Each 
annulus has 1 kpc width. Right: local definitions with 3" grid. The nucleus, bar, arm and inter-arm are represented 
by pink, green, red and blue, respectively. 
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Fig. 9. — : Histograms of (left) surface density Sh 2 , (middle) area averaged star formation rate Esfr, and (right) SFE. 
The insets are overall view of the histograms. 
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Fig. 10.—: Sh 2 vs ^SFr plot (Kcnnicutt-Schmidt plot). Left: Eh 2 vs Esfr at the resolution of 6" (~ 500 pc) (color) 
at the scale of Kcnnicutt's Figure 9 (1998). Right: the same as the left panel, but at a larger scale. 



